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The fusion process for heavy nuclei has been extensively studied within a gener- 
alized formulation of the diabatic approach to collective nuclear motion [1]. It was shown 
there that coupled equations of motion decouple during the approach if appropriate di- 
abatic channels are used. Thus the mean value and the variance of the fusion barrier 
are completely determined from the ensamble of diabatic barriers being defined by initial 
correlations. The mean shift AB with respect to the adiabatic barrier and the variance 
o} have been calculated from the splitting of the diabatic single particle levels at the 
barrier and the initial occupation probabilities given by the pairing wavefunctions of the 
separated nuclei. 

The microscopic model applied in [1] reproduces both the general (smooth) depen- 
dence of the mean shift AB on the mean fusibility parameter z,, and the strong opposite 
isotopic dependence due to shell structure present in experimental data [2,3]. The accu- 
racy of the theoretical predictions decreases, however, with the increasing fusibility. This 
is a signature that the existing theoretical model does not contain a physical mechanism 
which becomes important for heavier systems and sets on for z,,2/0.75. 


"Work partially supported by Grant IF/II/R . 


242 P. Rozmej, W. Nórenberg 


The principal results of the paper [1] ure presented in figure 1. 
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Fig. 1. Experimental (squares) and theoretical barrier shifts for 23 fusion systems as 
functions of the mean fusibility. The full dots give the values for the pure diabatic shift. + 
The extra push obtained from the one-body dissipation model [4] and the surface-friction 
model [5] are shown hy the dashed and dash-dotted lines, respectively. 


The model used in [1] can be considered as the diabatic limit, Barrier shifts and 
fluctuations are calculated under the assumption that the single-particle motion during 
the approach phase of the collision is purely diabatic. For heavier systems, as an overlap 
of the colliding nuclei at the barrier distance increases with fusibility, one can expect 
the increasing influence of two-body collisions. These collisions produce dissipation that 
generally broadens the barrier distribution. The mean shift is also increased because the 
additional energy ('diabatic extra-push') is required by the system to reach the fusion 
configuration. Two-body collisions are incorporated in DDD theory (Dissipative Diabatic 
Dynamics) [6]. The present paper aims to formulate DDD in a way convenient to describe 
the fusion process of heavy nuclei. As there is no single fusion barrier and contrary there is 
a definite distribution of fusion barriera (with the mean value and the variance calculated 
in [1]) we introduce reaction channels occuring with a definite probability. For simplicity 
we limit our considerations to only one collective degree of freedom q (in our case it is the 
distance between centers of colliding nuclei). Then the equation of the collective motion 
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for a particular channel v reads as 


ait? (2) e» 9. (a) 


where B is mass parameter and F,(q,t) is the force for the given channel. At the time 
t = 0 the force is given by the derivative of the potential that is the aum of adiabatic part 
and diabatic shift d 

Flat = 0) = -F [Vo Ate (2) 


The adiabatic potential V*4 determining the equilibrium force is taken from the frozen 
density Hartree-Fock (FDHF) calculations [1]. This paper supply also the distribution 
of the diabatic shift of the fusion barrier AV“*>. Choosing the different values of the’ 
diabatic shift one can define a particular reaction channel that is populated with the 
determined probability. Figure 2 shows the example of initial distribution of the fusion 
barriers obtained in static approach [1]. The mean value represent the sum of the adiabatic 
(FDHF) barrier and the mean value of the diabatic shift. 
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Fig. 2. Initial distribution of fusion barriers for the 199 Mo+™°Pd system. The thick solid 
line displays the mean barrier and the dashed and dotted lines correspond to barriers 
different from the mean value by one and 2 standard deviations, respectively. 
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Due to two-body collisions the force F,(q, t) fluctuates and is driven to its thermal 
equilibrium value. We may describe the stochastic behaviour of the force by a distribu- 
tion function which satisfies a Fokker-Planck equation. Alternatively we introduce the 
stochastic variable F,(g,t) which satisfies the Langevin-type equation 


ôF (g, t) 1 
"ee [EaD - FQ] +h), (3) 
where 7,, denotes local equilibration time. 


This time has been estimated by Bertsch [7] within Fermi-gas model as 


2.10775. MeV 
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and depends on the excitation energy e" per particle (which is related to the temperature). 
F(q) stands for the adiabatic force. The fluctuating force F#®*(q, t) is an irregular function 
in time with zero mean and a white spectrum (5 — functions in time) [8]. It can be related 
to the diffusion coefficient D of the corresponding Fokker-Planck equation [9]. For the 
actual numerical calculations we discretize equation (3) following [9] and obtaining 


Ers d Ee -— [£47 | r+ VD vr. (5) 


Here r is a time step and w(t) is a gaussian distributed random number, such that (w) = 0 
and (uw) = 2. 
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is the diffusion coefficient which is related to the mean value of the variance of the force 
2$. The last quantity is evaluated in the right hand side of (6) microscopically according 
to [8] and presented in figure 3. It appears to be closely related to stiffness parameter Ce 
(tensor for more collective coordinates) introduced within dissipative diabatic dynamics 
[10,11] via equation oi = T - C(T). Alternatively stiffness tensor is expressed as 


eres.) o 


The e, are the diabatic single-particle energies, ys is the chemical potential and 71, (q, 4, T) 
are equilibrium occupation numbers according to Fermi distribution : 


fi, = fı E enz Eil (8) 


Figure 3 displays stiffness parameter (sum of the proton and neutron contributions) for ' 
the system 1°°Mo+™°Pd calculated from diabatic single-particle states obtained in two- 
centre shel model which parameters had been adjusted to reproduce initial Hartree-Fock 
single-particle energies [1]. 

Numerical calculations of the fusion probability are not finished yet for the dy- 
namical model. Preliminary results (which estimate the mean dynamical shift of the 
fusion barrier) idicate that the dynamics of the process can indeed be responsible for the 
previous discrepancy between the theoretical and experimental results. 
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Fig. 3. The stiffness parameter for 189 Mo--!1? Pd system at different temperatures. Dotted, 
dashed, dash-dotted and solid lines correspond to temperatures T = 0.2, 0.6, 1.0 and 1.4 
MeV respectively. 
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